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SUMMARY 


A general analysis is presented which allows the rolling motions of 
an aircraft using a displacement— response, flicker— type automatic pilot 
to be determined. It is shown that the inertia, damping, and control 
characteristics of the aircraft in roll and the lag time of the automatic 
system are sufficient to define these rolling oscillations, and charts 
are presented by which the amplitude and period of the resultant steady- 
state oscillations may be found. Because of the inherent residual oscil- 
lations, this system is not considered ideal for many stabilization prob- 
lems; however, this flickei^-type system may offer a simple and economical 
solution to those applications where steady— state oscillations may not 
be objectional. Current trends in pilotless-aircraft designs indicate 
that the topic system can possibly provide roll stabilization with small 
amplitude residual oscillations. The analysis permits the definition of 
stabilization boundaries which will reveal whether or not a specific 
installation may prove satisfactory. A method for finding the transient 
conditions leading to steady state and the case of out— of— trim moments 
producing roll are also considered. Close agreement existed between the 
theoretical results and roll— simulator tests. 


INTRODUCTION 


Since the beginning of investigations into the field of pilotless- 
aircraft research, one of the most important and difficult of the funda- 
mental problems demanding solution has been that of roll stabilization 
by means of automatic control. Although considerable work has been done 
in the analyses of various automatic pilots, these efforts usually con- 
cerned the application of various systems to particular aircraft. This 
paper, however, represents an effort to study a single type of automatic 
pilot and its adaptability in terms of familiar aircraft parameters. 

Of the numerous possible methods for obtaining automatic roll stabi- 
lization, the one which has been selected for discussion here is the 
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displacement— response, flicker— type system. In this system the sense of 
the control moment is dependent solely on the sense of the angular 
displacement from a zero reference. The control moment is considered to 
he constant in either one direction or the other at all times. It is 
recognized that such a system may not he ideal for many stabilization 
problems; however, the flicker-type system may offer a simple and eco- 
nomic solution to those applications where residual oscillations of small 
magnitude are not objectionable. 

It is the purpose of this paper to determine how the aircraft and 
automatic-pilot parameters affect the rolling motions of an aircraft 
employing the topic system. This analysis includes the determination of 
the limits of the parameters involved in defining practical stabilization 
boundaries for the system. In addition, charts are presented which may 
be used independently from the mathematical analysis to determine the 
amplitude and frequency of the rolling oscillations of any given air 
craft. Eesults of tests of a subject automatic pilot and roll simulator 
are included to show the close agreement of this theoretical approach 
with experimental results. 


SYMBOIS 
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At 


angle of bank, radians 

angular rolling velocity, radians per second (d0/dt) 

p 

moment of inertia about longitudinal axis of aircraft, slug-feet 


rolling moment, foot-pounds 

rate of change of rolling moment with angular rolling velocity, 
footpounds per radian per second (dL/Bp) 

combined differential deflection of ailerons, radians 

rate of change of rolling moment with aileron deflection, foot- 
pounds per radian (SL/SS) 

control moment, foot-pounds 


out— of— trim moment producing roll, foot-pounds 


ratio of out— of— trim moment to control moment 


( | l o/ sl 6 | ) 


time^ seconds 

time increments measured from translated time origins ^ seconds 
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time lag between signal reversal and instantaneous control 
deflection, seconds 

damping-to— inertia ratio, per second ^|Lp/l x Q 

automatic stabilization parameter, dimensionless (st) 

maximum angular rolling velocity when out — of— trim and control 
moments are additive, radians per second (see equation (4)) 

maximum angular rolling velocity when out — of— trim and control 
moments are opposed, radians per second (see equation (22)) 

fractions of p^y 

fractions of p' mwy 

fraction of initial rolling velocity of a cycle existing at end 
of first half cycle 

fraction of initial rolling velocity of a cycle existing at end 
of complete cycle 

amplitude of steady— state rolling oscillations (one— half total 
displacement), radians 

displacement of mean line of steady— state rolling oscillations 
due to out^of— trim conditions, radians 

dimensionless amplitude factor 



period of the steady— state rolling oscillations, seconds 
natural logarithm 

base of natural logarithms (2. 71^3 ) 

basic change of angle of bank induced by control action, radians 
(see equation (13)) 

wing span of aircraft, feet 

wing area, square feet 

mass density of air, slugs per cubic foot 

velocity of aircraft, feet per second 

dynamic pressure, pounds per square foot (pV 2 /2^ 
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rolling— moment coefficient (L/qSb) 
nondimens ional damping— in— roll coefficient 



'*6 


rate of change of rolling-moment coefficient with aileron deflec- 
tion (SC j/d&) 


Numerical subscripts refer to time limits. (See fig. 1.) 


DESCRIPTION OF THE AUTOMATIC PILOT 


The displacement— response, flicker^type automatic-pilot design con- 
sidered in this paper as a roll-stahilization system consists primarily 
of one displacement gyroscope, a power supply, a servomotor, and an aero- 
dynamic control surface. In general, the analysis will apply to any 
automatic pilot which causes the application of a constant control moment 
for a displacement deviation. Such a system represents one of the simpler 
and more economic nonlinear servomechanisms used for automatic control. 

A schematic diagram of a typical system of this type is shown in 
figure 2. In operation, a rider arm is attached to the outer gimbal of 
the displacement, gyroscope which provides the constant zero space refer— 
ence. A commutator drum is attached to the mounting case of the gyroscope 
and provides the aircraft displacement reference. Thus, a deviation in 
the angle of bank is detected by contact of the rider arm upon the commu- 
tator drum. The commutator drum is divided into two segments, and contact 
of the rider arm with either of these segments completes a power supply 
circuit which energizes the servomotor. Energization of the servomotor 
causes an instantaneous deflection of the aerodynamic control and main- 
tains this deflection until contact is made with the other commutator 
segment . 

The intelligence of the system is such that the sense of the control 
moment is directed to restore the aircraft toward the zero space refers 
ence. Hence, as the aircraft rolls to the right, the rider arm is on the 
segment which completes the circuit calling for left control. In this 
analysis, controlled motions involving bank angles greater than l80 are 
not considered. 


METHOD OF ANALYSIS 


The analysis of the rolling motions of an aircraft controlled in 
roll by a displacement-response, flickei^-type automatic pilot can be 
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briefly summarized in three steps. First, a description of the angle— of— 
hank variation during a complete cycle of a rolling oscillation is written 
in terms of the general aircraft parameters which define the motion. 
Secondly, from the general expressions of the rolling motion, the rolling 
velocities at the beginning and end of the cycle are related; from this, 
the analysis shows whether the aircraft is performing a series of steady 
oscillations or whether the oscillations in the transient state are 
growing or damping to the steady— state conditions. The third step is 
the determination of the general expressions which define the amplitude 
and period of the rolling oscillations. 

Of primary importance in this analysis is the basic response of the 
aircraft to an instantaneously applied control moment. It is assumed in 
this analysis that motions in pitch or yaw do not affect the rolling 
behavior. Therefore, the basic response of the aircraft in roll to an 
abrupt application of control is given by the rolling-moment equation for 
one degree of freedom. 

Two mathematical principles are employed in representing the motion. 
The first of these principles allows the basic response of the aircraft 
to be redefined from different origins on the time scale; this feature 
is termed the principle of redefinition through time translation. The 
second principle is that of superposition of individual effects to detei^ 
mine resultant motions. The use of these principles becomes clear as 
actual applications are made. 


Basic equations of motion.— The basic differential equation of 
motion in roll in this analysis is the rolling-moment equation of 
one degree of freedom: 


The term L Q is included to account for any continuous out— of— trim 
moment of the aircraft that may be producing roll. The solution of this 
linear differential equation yields the relationships 


DERIVATION OF EQUATIONS 



( 1 ) 
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p . S1 6 4 L ° [ e (Lp/!x)t _ J + (3) 

Lp 


where 0(0) and p(0) are, respectively, the angle of bank and rolling 
velocity at the time origin considered. 

Equations (2) and (3) are modified in form for use in this paper. 
Since, in this analysis, no particular time is considered zero time 
because of the use of time translation, the notation in equations (2) 
and (3) is changed from t to At. In this case At indicates time 
increments from a defined origin, which is where At is considered 
equal to zero. Hence, equations will represent the variations following 
the translated time origins. A second change is made through the substi- 
tut ion of 



where a is termed the damp ing— to— inertia ratio and the vertical "bars 
indicate the absolute magnitude. A final modification is the expression 
of p(0) as a fraction of the maximum rolling velocity. From equa- 
tion ( 3 ) 


Pmax 


+ L o 


and 


P(°) = C Pmax 


(M 


(5) 


where C is the fraction of p^ at the time origin considered. 

Making these modifications, the basic equations of motion are 

0 = 5L & ^ |- (1 _ c)e ~aAt + a At _ (1 - C)] + 0(0) (6) 

ip 2 


&L& + L ° [l - (1 - C)e' aAt ] 

Lp 


P 


( 7 ) 
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First half— cycle.— In this analysis the out— of— trim moment L Q is 

assumed to he a constant positive moment causing roll to the right. At 
point 0, figure 1, the angle of hank is zero and the rate of roll is 
defined as Cn Pmflv With these initial conditions the angle— of— hank 

variation between points 0 and 1 is giver hy equation (6) as 


001 = S ~— g — ^ [(! - C 0) e_aAt + a At - (1 - C 0 )] (8) 


At zero angle of hank (point 0) the automatic pilot signals for a 
reversal of the controls. Because of the inherent lag in the operation 
of the physical system, the abrupt deflection of the control occurs at a 
finite time after the signal is given. As the aircraft is assumed to he 
rolling at this point, it reaches a definite angle of hank and rolling 
velocity at the end of this lag period (point 1, fig. l) . This lag time 
is considered constant for any system and is defined as t. Substituting 
this value into equation (8) gives the angle of hank at point 1 as 


01 = — p " [C 1 - C o) e_ &T + ax - (1 - C 0 )] (9) 

V 


The rolling velocity at point 1 is found to he 


Pi ^l^max 

where 

Ci = [l - (1 - Co)e~ aT ] (10) 

Tlie factor aT is of prime importance in this analysis. The 
symbol K will he termed the automatic stabilization parameter and is 
defined hy 



K = aT = 


T 


(ID 
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To completely define the result of the control reversal at point 1, 
it is necessary to define the angle — of— bank variation that would have 
existed had the control not reversed. The definition of this curve, 
shown dashed in figure 1, is accomplished by time translation. . If time 
increments are considered to be measured from point 1, the basic response 
equation of the aircraft may be used with a change in the initial coi>- 
stants, the new constants being the angle of bank and rolling velocity 
that were determined for point 1. In this way the rolling motion that 
would have occurred if the control had not reversed is redefined in terms 
of time increments from the time of the control reversal, which is con- 
sidered a new time origin. This is written as 



— - + _— lT(l - C 1 )e -aAt + a At - (1 - C-l)] + 

h > 2 


( 12 ) 


The effect of the reversal of the control in producing roll in the 
opposite direction is determined by superposition. The algebraic summa- 
tion of the motion that would have existed had the control not reversed 
with the incremental changes in angle of bank caused by the control 
reversal gives the resultant motion of the aircraft following the control 
reversal. 

The incremental changes introduced by the control reversal may be 
expressed by the basic response equation of the aircraft. The basic 
rolling moments produced are due to control deflections from the neutral 
to the maximum positions. However, at point 1, the control reversal is 
from a maximum deflection in one direction to the equal value in the 
other direction. Consequently, the incremental changes in the angle-of- 
bank variation introduced by the control action at point 1 are twice the 
magnitude given by the original basic response equation. The constants 
of integration for the control effect are determined by conditions of 
zero angle of bank and zero rolling velocity . This control effect is 
written as 


20 c = 2 — (e -aAt + a At - l) (13) 

h 

The resultant rolling motion between points 1 and 2, 0 12 , is found by 

the superposition of the motions expressed in equations (12) and ( 13 ) . 
This gives 
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01 2 = 5L& ! — ° ^ [(I - C 1 )e- aAt + a At - (l - C x )] 

h 

+ 0 , - 2 — 3_ (e"®** + a At - l) (lM 

r± t 2 X 


It is now defined that after a time t 12 seconds from the time of 
the control reversal the aircraft has reached a rate of roll of — 
or has reached the fraction C* of the rolling velocity existing at the 
"beginning of the cycle. The negative sign is used since the direction of 
the rolling velocity is opposite that defined at the signal reversal at 
point 0. This substitution into the rolling-velocity expression obtained 
from differentiation of equation (l4) gives 


-C’Cr 




'SLj= 







(1 - Cl )e^ t12 ] - 2 ^ (l - e^ 12 ) 

(15) 


Before solving equation (15) for t 12 it is convenient to make the 
following substitution: 


e = 


Lo 

6L 5 


(16) 


where e is termed the out-of-trim ratio. Substituting equation (l6) 
into equation ( 15 ) and using equation (10) for gives t^ 2 as 


J 12 


-a 


log y 


(17) 


where 



(18) 
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Furthermore , let t^_2 also he the time at which 0 = 0. (See fig. 1.) 
With 0=0 and At = t^ in equation (l4), and with equation ( 17 ) 
for t 12 , equation (l6) for L 0 , and equation (10) for C^, the fol- 
lowing form is reduced: 


(1 + e) [(l - C 0 )e -K ]Y - 2Y + (l - e) log Y + 2 + (l + e)(K - 1 + C 0 ) = 0 

(19) 


where Y is defined by expression (l8). 

Equation (19) is the conditional equation for the first half-cycle. 
This expression permits the determination of the rolling velocity at the 
end of the first half-cycle C l C n p Tn „ T when that at the start of the 

cycle CoPmax is known J these velocities are functions of the automatic 
stabilization parameter K and the oub-of— trim ratio e which must he 
known or estimated for the aircraft. 

Second half— cycl e.- At point 2 (fig. l) it is noted that the control 
moment is to the left and hence would require a negative sign. The sense 
of the out— of— trim moment is, of course, unchanged. In this case the 
result of this opposition of rolling moments is a change in the expres- 
sion that is given to the basic response equations of the aircraft. 

These equations are rewritten as 


where the maximum rolling velocity when the control and out— of— trim 
moments are opposed takes the form 


0 = ~ §L5 + L ° I x [(1 - c)e -aAt 4 a At - (1 - c)] + 0(0) (20) 



( 21 ) 


-8L S + L o 


( 22 ) 


max 



and c is defined as the fraction of 
considered. 


at the time origin 


NACA RM No. L8D19 


11 


The angle— of— hank variation between points 2 and 3 is IlOW required. 
At point 2 } 0 = 0 and the rolling velocity fraction c 2 for use in 

equations (20) and (2l) takes the form 

c 2 = c,c o(rfr) (23) 

when equation (l 6) is used for L Q . The variation 0 2 ^ is now written 
from equation (20) as 


023 - 


""Si'S + L o 

^7“ 


(1 - c 2 )e _aAt + a At - (1 - c 2 ) 


(2k) 


The steps used in the consideration of the first half-cycle are 
followed for the second half— cycle with the following important relatioi^- 
ships resulting: 


0, = 


-SL 5 + L o 


(1 - Cgje^ + K - 


V 



(25) 


p 3 C 3 P max 

c 3 = [l - (1 - c 2 ) e ^ 

The motion 0 2 ^ which would have resulted had the control not reversed 
at point 3 is given as 

02 3 4 ' ~ 61 ' 5 y - ^[i 1 - <= 3 ) e_aAt + a at - f 1 - c 3i 

V 


(26) 

(27) 


+ 03 (28) 
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and by superposition the resultant rolling motion 0^4 following the 
control reversal takes the form 

(*31, - ^3 L ° ^[(1 - <=3) 8_aAt + a At - (1 - o 3)] 

+ + 2 — + a At - l) (29) 

h 2 

Expressing the rolling ‘velocity after a time t^i^ seconds from the 
second control, reversal as C M CoP m R -r j or as the fraction C M of the 
rolling velocity at the beginning of the cycle ^ allows t^ to he 
written as 


t 34 - 


— log 


where 


X = 



(30) 


(31) 


Finally j the conditional equation for the second half-cycle is reduced 
to the form 


(1 - e) (l — X) 


1 - C'Cq 


1 + e 


.-K 


- (1 - e) < — log X + K 


(e" K - l) 


1 _ c'Cofj-^ 


+ 2(X - log X - 1) =0 (32) 


Equation (32) allows the determination of the rolling velocity at the 

end of the cycle C"C„p when the rolling velocity at the start of 

u max 
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the second half-cycle c ' c oPmax> the automatic stabilization parameter K, 
and the out— of— trim ratio e are known. 

Amplitude equation.— The amplitude of the oscillation is defined as 
one-half the total displacement. The total displacement is found from 
the maximum angles of bank for each half-cycle. Since the maximum angle 
of bank: occurs at zero rolling velocity, the substitution of the time to 
reach zero rolling velocity into the angle— of— bank equation for the 
desired portion will give the required maximum angle. These results 
are written as 


(^isW ” B ' 


(1 + e ) (k + C 0 ) + (1 - e) log 


1 - e 


2 - (1 + e)(l - C^e-K 


(33) 


= “ B ( Kfl “ e) + C ‘°0( 1 + 


+ (1 + e) log < - 


1 + 6 


2 - (1 - e )e“^ 


- c '°o(yr 


( 34 ) 


where B is termed the amplitude factor and is defined as 




V 


The amplitude A of the oscillation is given as 


(35) 


A = 2^12)^ - (^naj 


(36) 


and is calculated by using equations ( 33 ) and ( 3 * 0 . 

Mean-line displacement.— Since the presence of the out— of— trim 

moment causes the oscillation to be unsymmetrical about the zero angle— 
of— bank position, the mean line of the oscillation is displaced from zero 
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in the direction of the out— of— trim moment. The displacement of the mean 
line is designated Aq and is given by the equation 


= If^max + (^)inax] 


(37) 


Period equation.— The period of the oscillation is defined as the 


time required to complete 1 cycle. From figure 1 it is evident that 
the period P may he expressed as 

P = 2T + ti r\ + t 


12 + 


(38) 


Substituting equation (17) for t 12 and equation (3O) for t^ allows 
equation (38) to he written as 


P = r [2 - i(lOg X + log Y)] 


(39) 


where X and Y are defined hy equations (31) and (l8), respectively. 

Angle— of— hank: equation.— The angle-of— bank time equation is written 
in four parts . Each part is computed separately, and time increments 
are from At = 0 for each part. The equation for the resultant angle- 
of— hank variation is 


0 = 0Q1 + 012 + 023 + 034 


(40) 


where the equations are (8) for 0 qi* (l4) for 0^2* (24) for 023* 

and (29) for 0^. 


PRESENTATION OF RESULTS 
General Remarks 


Prior to the application of the derived equations to the detailed 
analysis, several observations of some of the actual physical properties 
of an aircraft-automatic pilot combination are presented. In this 
instance, as in nearly every other system, it is the actual physical 
characteristics that define the limitations of the automatic stabilization 
system. 
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In any real automatic-pilot system of the type under discussion, 
there is in its operation a definite time lag — the time between the 
sensing of an error and the application of corrective control. Con- 
sidering the statement above, it is readily possible to include all 
forms of lag, whether they be considered mechanical, electrical, aero— 
dynamic, or otherwise. Because of this lag the aircraft, when once 
disturbed, must oscillate, these oscillations never damping completely. 

It Is only in the theoretical case of zero lag, which represents a limit 
of possibilities for this system rather than a possibility itself, that 
the induced oscillations following a disturbance will damp to zero 
amplitude . 

In general, it may be stated that within certain ranges of various 
parameters the automatically controlled aircraft will, upon being dis- 
turbed, tend to oscillate at a definite amplitude and frequency. This 
condition is defined as the steady— state condition, and at this steady- 
state condition the aircraft will never exceed a definite rate of roll 
which is real and determinable. If the disturbance is such as to cause 
the rate of roll at the time of signal reversal to be greater than the 
steady— state value, the transient oscillations which follow will damp 
to the steady— state conditions; and if the disturbance is such as to 
cause the rate of roll to be less than the steady— state value, the 
resulting transient oscillations will build up to the steady— state values. 
This steady— state condition is a real limit which theoretically is only 
approached; however, for all practical purposes* it is usually very 
closely approached in a relatively short time after the disappearing 
disturbance is encountered. 


Transient and Steady State 

The conditional equations, defined by expressions (19) and (32), 
present the relationships which will describe the actions of the aircraft 
under any prescribed condition. The prescribed condition necessarily 
involves the automatic stabilization parameter K and the out— of— trim 
factor e; for a given aircraft these parameters must either be known 
or estimated. The other factor of the prescribed condition is the 
fraction Cq of maximum rolling velocity at the time of a signal 
reversal. 

The conditional equation for the first half— cycle will show the 
fraction C f of the initial rolling velocity considered c oPmax 

exists at the time of the second signal reversal, provided e and K 
are given. The use of the second conditional equation presupposes 
that C f has been determined by use of the first conditional equation, 
and it reveals the rolling velocity which exists at the end of the 
cycle C"C oP[1M . 
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Transient state.- The two conditional equations may "be used to show 
the transient conditions occurring as the aircraft approaches the steady 
state. An evaluation of the two conditional equations which reveals C" 
to he less than unity shows that the rolling velocities of the aircraft 
are decreasing, and hence, the aircraft motions are damping to the 
steady— state condition. The increase to the steady— state condition is 
shown hy values of C" greater than one. 

As an example, the analysis is applied to the determination of 
transient conditions with the assumption of perfect trim (e = 0) . In 
this case the rolling oscillations are symmetrical about zero angle of 
hank. For a value of K = 0.5, the respective values of C' and C" 
were computed for various values of C 0 hy using the two conditional 
equations. These results are presented as figure 3 with C" plotted 
against C 0 . Although figure 3 is a particular case, the shape of this 
curve is typical of those obtained when any values of € and K are 
used. 


As an illustration of the use of this type of curve, let a value of 
Cq = 0.2 he considered. From figure 3, "the rolling velocity at the end 

of that cycle C w C n p Tt|aT is shown to he (3*55) (O^P^y = 

Considering C 0 = 0.71 for the next cycle, the rolling velocity at the 

end of this second cycle is given from figure 3 again as 
(l.05)(0.7l)p nax; = (0.7 1 <-55)p maz - Continuing this method would reveal 

that the condition of C" = 1.0 and C 0 = 0.75 is approached; this is 


the approached steady— state condition since all of the cycles become 
approximately identical. Assuming a value of Cq greater than the 
steady— state value will also show the tendency to approach the steady- 
state condition. Thus, when C" < 1, the damping character of the 
motion is shown, and the increasing action is revealed when C" > 1. 


Steady state.— At the steady-state condition, the value of C" must 
he unity since this case assures the equivalence of every cycle. Since C 
is defined when Cq is defined (equation (19)), the steady— state condi- 
tions may he expressed solely as functions of 6, K, and Cq. Equa- 
tions (19) and (32) have been solved for steady— state conditions and the 
results presented in figure U as K against Cq with e as a parameter. 
Thus, for given values of e and K, the value of Cq at steady-state 
is shown hy figure 4. 


Stabilization Boundaries 

In this analysis motions involving angles of hank greater than 180 
are not considered. This consideration allows limiting values of various 
parameters to he defined* which will permit this maximum angle-of-bank 
condition. Indeed, limiting values can be defined for any maximum angle 


NACA EM No. L8D19 


17 


of 138X115:. These upper limits of the various factors are termed stabiliza- 
tion boundaries. These boundaries are presented to allow a rapid evalua- 
tion of whether or not a given aircraft employing this automatic pilot 
may possibly give satisfactory amplitude characteristics. 

Since the effects of aut>~of— trim conditions in producing roll are 
included, the maximum angle of bank will be in the direction of the out— 
of— trim moment (to the right in this analysis). This maximum angle of 
bank has been defined in equation ( 33 ) as ( 012 ) m7 • Substituting 

jr radians ( 180 °) into this expression will give the relationship that 
must exist between the parameters to insure this condition. This expres- 
sion is 


rt = B 


(1 + e )(K + C 0 ) + (1 - e) log 


1 - e 


2 - (1 + e)(l - C 0 )e- K _ 


(hi) 


It is noted that, at steady state, the bracketed term in equation (kl) 
is defined when € and K are defined. Figure 5 is presented as the 
solution of equation (4l) and is expressed as a function of the amplitude 
factor B and the automatic stabilization parameter K for various 
out— of— trim ratios e. Also shown on figure 5 is the boundary where the 
limit maximum angle of bank is n /6 radians ( 30 °). 

If the values of K and B for an aircraft are plotted as a point 
on this graph and the point is below the boundary, the aircraft will 
stabilize in roll at an amplitude less than the value shown for the 
boundary. These curves indicate clearly that the magnitude of these 
factors is rather limited. 


Amplitude Equations 

The amplitude of the rolling oscillations, defined by equation ( 36 ), 
is considered herein only for the steady state. Since values of Co* 

C f , and C" are defined at steady state when e and K are known, 
the amplitude A may be considered as a function of three parameters, 

€, K, and B. Since the steady— state amplitude A varies directly 
with the amplitude factor B, it is possible to present A/B against 
the automatic stabilization parameter K for any given value of the 
ou t— of — tr im rat i o € . 

Figure 6 shows the A/B variation with K for the case of perfect 
trim ( e = 0) . (Value of A computed from this figure is in degrees.) 

The consideration of oui>-of— trim conditions led to calculations of 
the steady— state amplitudes for several values of e over the range 
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of K from 0 to 4. These calculations show that at any value of e 
the value of A/B was always within 6 percent of the value of A/B 
at e = 0 for the same value of K. In other words , the amplitude of 
the steady— state rolling oscillations (one— half of the total displace- 
ment) is practically unaffected "by out— of— trim moments. 

In understanding this action, it is necessary to remember that K 
and B are constant and only 6 is changed. With the introduction 
of e, it is seen that the unsymmetrical character of the oscillations 
appears, the presence of the out— of— trim moment causing the control 
action to be slow in stopping roll in one direction and very quick in 
stopping roll in the other direction. This effect, however, results in 
very slight changes in the total amplitude displacement since the out-of- 
trim condition causes a greater maximum angle of bank in one direction 
but a lesser maximum value of bank in the other direction. Consequently, 
figure 6 may also be used satisfactorily in estimating the amplitude of 
the steady— state roll oscillations under out— of— trim conditions. 

Mean-line displacement .— The unsymmetrical nature of the rolling 

motions under ou"t>-of— trim conditions may be expressed as the displacement 
of the mean line of the oscillations from the zero reference. Equa- 
tion (37) represents this displacement which is rewritten here as 

A 0 = | [f 12(6, K) + f 3 k(e, K)] (42) 

where f^fe, *0 and f^(e, K) are the "bracketed terms of equations (33) 
and ( 34)3 respectively ^ "both of which are completely defined hy c and K 
under the steady-state conditions which are under consideration here. 

The solution of equation (42) for the steady state is presented in 
figure 7 . The ratio of the mean-line displacement A Q to the amplitude 
factor B is expressed as a function of e and K. In order to deter- 
mine the extent that the mean line of the rolling oscillations has been 
displaced ^ the value of the ordinate of figure 7* at the values of K 
and € in quest ion, multiplied by the amplitude factor B will give 
the value of the displacement A Q in degrees. 


Period Equations 

The determination of the period of the rolling oscillations is 
accomplished by use of equation (39). Figure 8 presents the graphical 
interpretation of equation ( 39 ) for steady-state conditions , and is a 
plot of the ratio of the period P to the lag time t against K for 
various values of e. Using the value of K for the aircraft 9 the 
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P/t value read from figure 8 multiplied by the lag time t will give 
the period of the steady— state oscillations. 

Although the ratio of the period to the lag time is presented, care 
should he taken in interpreting the significance of this ratio. The 
period does not vary directly with the lag time since K is also a 
function of t. For example, the doubling of t does not double the 
period since the value of P/t is changed due to the resultant doubling 
of K. 

The dotted portions of the curves of figure 8 are extrapolations. 

This was a result of the determination of the natural logarithm of the 
factor Y in equation (39) which involves a very careful consideration 
of the extent to which this decimal fraction is evaluated. A more exact 
solution was not considered Justified, particularly under the more extreme 
out— of— trim conditions. 


APPLICATION OF RESULTS 


It is evident that the results of this analysis can be used to 
determine the steady— state amplitude and period without actually substi- 
tuting into equations. This procedure may be broken down into five steps. 

(1) Determine the factors K, e, and B. 

(2) Plot K and B as a point on figure 5. 

( 3 ) Calculate the amplitude from figure 6. 

(4) Calculate the displacement of the mean line from figure 7. 

(5) Calculate the period from figure 8. 

Step (2) will indicate immediately whether or not the conditions might 
possibly satisfy some specific amplitude requirement. The three latter 
steps will give the actual values at steady state. If the actual bank 
variation is desired, equation (40) should be used. 

In an effort to indicate whether or not this automatic roll-stabili- 
zation system might be of use in current problems, table I was prepared. 
The aircraft shown are not identical with any particular pilotless — 
aircraft designs, although current trends are revealed in the magnitude of 
the characteristics presented. (The silhouette of aircraft 5 is approxi- 
mately one— half the relative size shown for the others.) In the case of 
aircraft 3 , the effect of doubling the lag time is shown to increase the 
amplitude and period by a factor of about 1.7 under the conditions con- 
sidered. Aircraft 1 and 3 show the possible importance of the amplitude 
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factor B. In “both cases the values of K are comparable, hut because 
of B the amplitude for aircraft 2 is extremely large , whereas air- 
craft 1 may be satisfactory for some applications. 

Although no special trends can be determined from this table , it 
does indicate that satisfactory roll stabilization with the topic system 
is a definite possibility in some cases. The method presented herein 
allows a rapid determination of the rolling characteristics and for a 
definite aircraft would be of considerable aid in determining the 
effects of variations of individual parameters. 

In the following paragraphs the present method will also be shown 
for one of the cases checked by roll— simulator tests. Attention will 
also be called to other additional factors which affect the application 
of this method. 


Boll-Simulator Tests 

In order to substantiate these theoretical results, experimental 
tests were run using the roll simulator developed by the Instrument 
Eesearch Division of the Langley Laboratory. This instrument is a 
single degree of freedom system which is controlled by an actual auto- 
matic pilot mounted on the moving table. The control torque— to— inertia 
ratio and aircraft damping— to— inertia ratio are electromechanically 
simulated. 

An Azon gyroscopic unit and an electrically operated servomotor 
were used as a typical displacement— response, flickei^type automatic 
pilot. In these tests the values of the torque— to— inertia ratio (SL^/Ix) 

and the damping— to— inertia ratio (Lp/L^ were varied. Oscillograph 

records of angle— of— bank and control position (servo-operation) against 
time were taken. Figure 9 shows the recorded results for two of the 
cases tested. 

It is noted that the experimental control motion is trapezoidal in 
shape. It is apparent , therefore, that the lag time cannot be taken as 
the time period between the signal reversal at zero angle of bank and 
the time the control reached full reverse deflection because to do so 
would be to neglect the effect of the control during the time it was 
reversing. Since the theoretical analysis is based on the instantaneous 
application of the maximum control torque, an equivalent square wave was 
substituted for the control motion. An effective lag t was computed 
as the time between signal reversal (at zero angle of bank) and the time 
of control reversal given by the equivalent square wave control motion. 
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Two illustrative cases which were tested are presented herein. 


Case 1 .— Calculations using the present method are given for the 
conditions employed in case 1 ; the experimental results are presented 
in f igure 9 (a) . 

Given: 


6L5 


32.0 


a 


h> 


= 4.0 


Calculated: 




2.0 


T = 0.025 Bee 

K = aT = 4.0(0.025) = 0.10 

A 

From figure 6, at K = 0.10 and e = 0, — = 8.0. Therefore, 

B 

A = 8 . 0(2.0) = 16.0° 

The average measured amplitude equals 15.7°. From figure 8 , at K = 0.10 
and e = 0 ^ 

p 

j = 21.2 


P = 21.2(0.025) = O. 53 O sec 

The average measured period equals 0.511 second. 

Case 2 .— Case 2 represents typical conditions for an aircraft similar 
to number 3 (table I) at a Mach number of 0.6. Figure 9(b) shows the 
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simulator results. For this case values of 


6Lg 


= 43.5, a = 9.jk, 


and t = 0.026 were used. These data give a calculated amplitude 
of 8.95°, whereas the average measured amplitude is 10.8°. The calculated 
period is O .355 second as compared with an average measured period of 
0.409 second. 


The differences encountered in the measured and calculated values 
are well within the limits of accuracy afforded the simulator at the 
present time. Additional error is "believed to have entered case 2 due 
to the fact that the employed torque-to-inertia and damping-to-inertia 
ratios are out of the linear portions of the simulator characteristics. 


Additional Factors 


In discussing the topic of automatic roll stabilization, there are 
several additional factors that should he presented in order that the 
applicability of the preceding analysis may he more clearly realized. 
These factors concern some of the included parameters, their significance 
and use, as well as some others not specifically included in this paper. 


Effects of velocity and air density .— Two primary factors involved 
in the analysis of the rolling motions of an aircraft are affected con- 
siderably by changes in velocity and air density; these parameters 
are Lp and Lg. Both of these derivatives are in terms of actual 

rolling-moment variation and hence are affected by velocity and air 
density. It can be shown that, for a constant damping-in-roll coeffi- 
cient C 1 , the parameter Lp varies directly as the product of the 


2V 

air density p and the velocity V. It may also be shown that the 
control parameter Lg varies directly with the dynamic pressure q 


when the coefficient Cj 6 is assumed constant. 


It is necessary, therefore, to emphasize that the analysis herein 
will pertain only to those conditions upon which the determination of Lp 
and Lg are based. The effects of changes in velocity and/or air 

density on the automatic roll-stabilization characteristics of an air- 
craft require that they be investigated as desired by particular cases. 

Effects of compressibility .- In the previous paragraphs it was 
mentioned that the coefficients ^ and C 2 g were considered constant 

2V 

in determining the effects of velocity and air density. However, the 


NACA RM Wo. L 8 D 19 


23 


constancy of these coefficients, particularly in the transonic speed 
range, is not assured due to the compressibility effects. If automatic 
roll stabilization is contemplated in or through this speed range, the 
stabilization characteristics should be based on whatever theoretical or 
experimental data are available concerning the coefficients. 

Significance of 8Lg.— In this analysis the actual moment produced 

by the control is expressed in aerodynamic notation 5Lg. However, the 

corrective control moments may not necessarily be aerodynamic; for 
example, the use of intermittent rocket Jet blasts has been suggested 
as a possible control device, particularly for use on extremely high- 
altitude pilotless aircraft or guided missiles. It is emphasized that 
the term 8L§ used in this paper refers to the actual moment used in 
controlling the aircraft, regardless of its nature; however, care should 
be exercised to insure a proper application to the problem when systems 
other than aerodynamic controls are used. 


CONCLUDING REMARKS 


A general analysis of the rolling motions of an aircraft using a 
displacement-response, flickej>-type automatic pilot has been presented. 
Because of the inherent residual oscillations, this system is not con- 
sidered ideal for many stabilizations problems; however, this flicker- 
type system may possibly offer a simple and economic solution to those 
applications where steady— state oscillations may not be objectionable. 
Current trends in pilotless-aircraft designs indicate that the topic 
system can possibly provide roll stabilization with small amplitude 
residual oscillations. Limits of the parameters involved have been 
presented as stabilization boundaries which allow a determination of 
whether or not a definite installation may possibly be satisfactory. 
Charts have been prepared for the determination of the amplitude and 
period of the resultant steady— state oscillations, and these theoretical 
considerations have shown close agreement with roll— simulator tests. A 
method for finding the conditions of the transient oscillations leading 
to steady state has been outlined. Consideration has also been given to 
the case where out— of— trim moments causing roll are present, and this 
effect is discussed. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I 




[ 




■ 

| 



| 

'< 

1 


1 

■ 


Aircraft 

i 

2 

3 

4 

5 

Weight, lb 

600 

150 

500 

1200 

5200 

Diameter, in. 

10 

8 

20 

17 

34 

Length/Dia . 

23 

15 

6.0 

12 

9.6 

Aspect Ratio 

1.5 

2.5 

4.6 

4.0 

5.1 

Mach Number 

1.7 

1.5 

1.5 

0.85 

0.65 

6 L*. ft -lb 

213 

347 

930 

270 

600 

9- 1 

L , ft- lb/ rad/ sec 

Jr 

50.5 

6.58 

420 

17.5 

3890 

I x , slug-ft 2 

2.1 

0.3 

7.8 

14.3 

1665 

a = L^/l x 

24.0 

21.9 

54.0 

1.18 

2.34 

X f sec 

0.025 

0.025 

3.025 

a 050 

0.025 

0.025 

ii 

p 

0.60 

0.55 

1.35 

2.70 

0.0295 

0.058 

B = («%)i x Ap 2 

0.18 

2.41 

.3.041 

ao4i 

13.0 

0.066 

Fig. 6 A/B 

43.0 

39.5 

, 90.5 

171 

2.5 

4.8 

Fig. 7 P/r 

9.3 

- ii® 

I 6.7 

5.4 

30.0 

28.0 

Amplitude, ieg 

7.75 

95 

3.71 

7.0 

32.5 

0.315 

Period, sec 

0.232 

0.24 

3.168 

0.27 
1 

0.75 

0.70 
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